Escape from noisy intermittent repellers 
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Intermittent or marginally-stable repellers are commonly associated with a power law decay in the 
survival fraction. We show here that the presence of weak additive noise alters the spectrum of the 
Perron - Frobenius operator significantly giving rise to exponential decays even in systems that are 
otherwise regular. 
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There exist a variety of physical situations where one 
deals with the escape of trajectories from repellers. For 
instance, nuclear physicists are interested in the escape 
of particles along fission channels. Similarly, in the study 
of transport coefficients in two-dimensional ballistic con- 
ductors, one has to worry about the trapping time and 
its relationship with the geometry of the system 0] . By 
and large, it is now accepted that hyperbolic (chaotic) 
dynamics leads to an exponential decay in the number 
of trapped particles while intermittency or marginal sta- 
bility results in power law decays. Thus, a rectangular 
billiard table with a hole in the wall shows a power law 
decrease in the survival fraction while a (hyperbolic) en- 
closure created by the intersection of three discs shows 
an exponential decay in the number of trapped particles. 
The reason for this difference is intuitively clear from 
the following argument: Consider that there are N par- 
ticles distributed uniformly in (Birkhoff) phase space ||] 
and there is a hole along the wall at qo of extent Ago- 
The (average) fraction of particles that escape at each 
bounce is identical for a chaotic system and proportional 
to Ago since the particles remain uniformly distributed 
with time. Thus the survival fraction decays exponen- 
tially. In a marginally stable system however, an initial 
uniform distribution does not remain uniform at each 
bounce since individual particles tend to stick around 
stable islands. A heuristic derivation of the power law 
decay can be found in jsj and we merely remark here 
that the decay exponent is often difficult to determine 
analytically and an interesting advancement in this di- 
rection has been achieved recently by Dalhqvist Q . 

We are interested here in a situation where marginal 
stability or intermittency is accompanied by weak addi- 
tive noise. Such a situation can arise for instance in an 
imperfectly fabricated ballistic conductor in the shape of 
a triangle or stadium where reflection is no longer spec- 
ular but has additive noise. Thus : 



Qn+l ^ fliqn,Pn) 
Pn+1 = f2{qn,P,i) + 6i 



(1) 



where q,p are the Birkhoff coordinates ||], /i, /2 are the 
bounce maps and ^„ is a random variable with (^„) = 
having a normalized distribution g{^) (normally taken to 
be a Gaussian with zero mean). The question that we 
pose is : does one expect to find a power law decay in 
such a situation ? The answer we believe is interesting 



and can significantly alter the way people look at signa- 
tures of low-dimensional chaos in various experimental 
situations where noise is inevitable and often desirable. 
First, however, we shall consider a 1-dimensional inter- 
mittent map and study its spectrum in the presence of 
weak noise. 

A trajectory in the presence of additive noise is gen- 
erated by the iteration 



(2) 



where /(x) is a map, f„ is a random variable as described 
above and xq € [a,b]. An initial density of trajectories 
x) evolves according to the Perron - Frobenius equation 



(£oo0)(a;) = / dy6{x~f(ymy) 



(3) 



in the unperturbed case. Thus, the eigenvalues and eigen- 
functions of Co determine the escape rate in an open sys- 
tems. More specifically, assuming that the spectrum is 
discreet, an initial density 4>{x) can be expanded as 



c) = y^^Cgipajx) 



(4) 



so that the fraction of particles that survive n iterates of 
the map is 

p^^^ ^ /g dx (£g o (j)) jx) ^ ^ dxipa^ix) 



_ g-n ln(l/Ao) 



/a dx (t>{x) 



as n 



(5) 



In the above, {Aq} are the eigenvalues corresponding to 
the eigenfunctions {Lpa{x)} and Aq is the leading eigen- 
value with the largest real part. The discreetness assump- 
tion however holds only when the dynamics is hyperbolic. 
In the presence of marginally stable cycles, the spectrum 
has a continuous part leading to a power law decay of 
correlations (in closed systems) or survival fraction (in 
open systems). 

The presence of additive noise results in a modified 
kernel whose formal expression is well known [|| : 



{Co(j)){x)^l dy C{x,y)(p{y) 



(6) 



where 



1 



C{x,y) = Sa{x - f{y)) 
S„{x) = / S{x-Og{0 di^g{x). 



(7) 

(8) 



As before, if the spectrum of £ is discreet, the decay 
is exponential and the leading eigenvalue determines the 
asymptotic decay rate. 

Note that we have so far steered clear of spill-over 
effects due to noise 0. The most commonly adopted 
technique is the use of periodic boundary conditions which 
avoids spill-over altogether. Alternately, one can work in 
the infinite domain (—00,00) so that natural boundary 
conditions may be employed. Yet another approach is to 
tailor the noise distribution so that the probability of the 
dynamical variable escaping from the interval is zero. We 
shall, in this paper, have occasion to use the second and 
third approaches depending on the problem and it must 
be noted that there are other approaches to the spill-over 
problem that may be more realistic in a given situation. 
Note that the spectrum and eigenfunctions of C can be 
sensitive to the choice of boundary conditions. 

With this background, we now introduce the intermit- 
tent map 



(9) 



'x{l + x'^) x<0 
f{x) ^ ^x{l+p {2xy) < x < 1/2 
2a; -1 a; > 1/2 



where s > and p > I. The particle is considered to 
escape if a;„+i < or Xn+i > 1. The map is defined in 
the infinite domain so that natural boundary conditions 
apply on the density. The intermittency here is due to 
the fact that /'(O) = 1 so that the fixed point a; = is 
marginally or neutrally stable. 

For an initial uniform distribution of particles in [0, 1], 
the fraction that survives one iterate is clearly the sum 
of the two intervals II and Ir for which < f{x) < 1. 
Similarly, the fraction that survives two iterates is the 
sum of the four intervals III, IlRt Irl, Irr for which 
< /^(a^) < 1- Generalization in this case (of binary 
symbolic dynamics) is simple : the fraction that survives 
n iterates is the sum of the 2" intervals for each of which 
< /"(a;) < 1. Each of these intervals contains a peri- 
odic point and the larger its (in)stability, the smaller the 
size of the interval. Thus 1^ : 



j{n} ^ ^ 

A 



(10) 



where g is a symbol sequence of length n consisting of 
L and R which denotes the order in which the left and 
right branches (with respect to x = 1/2) of the map are 
visited, Oq is a constant and A^ — -^f'^{x)\x£i^ is the 
stability of the periodic point. The survival fraction can 
thus be expressed as [0,p| : 



{»} 



p r—1 



1 

27ri 



dz 



log C-'(^) dz 



(11) 



where ( ^{z) = Y[p (1 — z'^p / \Ap\) is the dynamical zeta 
function and 7 is a (small) negatively oriented contour 
around the origin. Dahlqvist [Q has recently shown that 
in the noiseless case, C~^(z) has a singularity of the type 
(1 - z)!/". It then follows 0] from eq. Q that the 
survival fraction, T{n) ^ for an initial uniform 

distribution of particles. 

The zeta function is also (approximately) related to 
the eigenvalues of C through the relation : 



Tr £" = 



^as-EEtt 

a p r—1 ' 



-2n 



(12) 



When the zeta function is analytic, its zeroes, {zk\ are 
isolated and related to A^ as A^ — l/zfc. On the other 
hand, when the system is intermittent and i!l~^{z) dis- 
plays a branch cut, the spectrum of C no longer re- 
mains discreet. Thus intermittency leads to an asymp- 
totic power law decay. For s = 0.7, p = 1.2, the initial 
decay is however exponential and this is ascribed to a pair 
of complex conjugate roots Q|. The power law behaviour 
emerges only after 600 iterations of the map. 

We now consider the map (eq.(^) with weak additive 
Gaussian noise 



5(0 



1 _^ 



(13) 



V2^2 

for s = 0.9, 0.7 and 0.5 and with a = 0.002 (see figure |) 




c 




200 400 600 ^800 1000 1200 1400 



FIG. 1. The survival fraction, T{n) (log scale) for s = 0.5, 
0.7 and 0.9 with a — 0.002. The straight lines (dashed) are 
the best fit exponential decays. The initial distribution in 
each case consists of 10^" points. 

Clearly there is a transition from a power law to an expo- 
nential decay in the presence of weak noise in each of the 
three cases for large n. Note, however, that the initial 
decay, though exponential, is at a significantly different 
rate and the slope of r(n) settles down to the asymptotic 
value gradually after a large number of iterations. 

Two inferences can be drawn from this transition from 
power-law to exponential behaviour. First, the presence 
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of weak noise makes the eigenvalues (of £, see eq. ^ ) 
discreet. Also, there are closely spaced eigenvalues with 
small differences in their real parts around the leading 
eigenvalue Aq which leads to the gradual change in slope 
of r(n). 

The discreetness of the spectrum follows from the fact 
that the noisy kernel is integrable and bounded. The 
corresponding Fredholm determinant ]lO| is thus entire 
whose zeroes (I/Aq) are isolated. Thus, even for very 
weak noise, the spectrum is discreet though the transition 
time may be too large for the final exponential decay to 
be observed experimentally. 

The closely spaced eigenvalues around Aq are possi- 
bly remnants of the continuous spectrum which exists for 
the noiseless case. In order to understand this better, we 
have evaluated the eigenvalues of £ by discreetizing the 
integral equation and diagonalizing the resulting matrix 
1^. Recall that a Fredholm integral equation exists as 
a limit of discreet sum fl0|] so that a matrix representa- 
tion is adequate so long as its order is large and spurious 
eigenvalues are eliminated. Figure ^ shows a plot of the 
eigenvalues of £ for s = 0.7 and a — 0.002 (see eq. (^)) 
obtained using a matrix of size 2500 x 2500. We have 
checked that the relevant eigenvalues on the positive real 
axis have converged to the fourth significant digit. 
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FIG. 2. The spectrum of £ for s = 0.7 and a = 0.002. 



Clearly, the closely spaced eigenvalues along the real line 
do not allow the survival-fraction to be dominated by the 
leading eigenvalue for small n. An order of magnitude 
evaluation of the transition time can be made by noting 
that T{n) ~ e-"''^(i/^«) {do + die'^^''^^"^^^'') where Ai 
is the next to leading eigenvalue and da — Ca J dx(j)a{x). 
Thus, for n» ntrans — l/ln(Ao/Ai), the leading eigen- 
value dominates. For the three values of s considered, 
ntrans — 98,105,119 for s = 0.5,0.7,0.9 respectively so 
that exponential decay sets in first for s = 0.5 as ob- 
served in fig. |l|. Thus the difference between the leading 
and the next-to-leading eigenvalue gives a good qualita- 
tive picture and fixes a lower bound for the transition 
time. Note also that for each of the three s-values con- 



sidered, the leading eigenvalue accurately reproduces the 
asymptotic decay of the survival-fraction (see Table |). 

We now turn our attention to the effect of noise in 
non-chaotic systems. Specifically, we shall consider trian- 
gular billiards which are non-chaotic though generically 
non-integrable. The only integrable examples are the 
(7r/3,7r/3,7r/3), (7r/2, 7r/3, 7r/6) and the (7r/2, 7r/4, 7r/4) 
triangles while all other rational triangles have at least 
one internal angle of the form rm:/n,m > 1 and 
are non-integrable. Their invariant surface, though 2- 
dimensional, is not a torus but topologically equivalent 
to a sphere with multiple holes. Also these systems are 
non-chaotic though irrational triangles are possibly er- 
godic and even display the weak mixing property [pT| . 
A linear stability analysis shows that the Jacobian ma- 
trix has unit eigenvalues and hence these billiards are 
marg inally stable Q. 

Consider such a triangular billiard of unit perimeter 
and let q,p denote the Birkhoff co-ordinates. Here q is 
measured along the boundary while p = sm(9) where 9 
is the angle between the ray and the inward normal at 
the boundary point q. Thus q G [0,1] and p € [—1,1]. 
In a typical experiment, one considers an initial uni- 
form distribution of particles (~ 10*) in this phase space 
evolving freely in between bounces and reflecting specu- 
larly from the walls. The particles are allowed to escape 
through a small opening at qo of extent Aqo{= 0.005). 
For both the integrable (7r/2,7r/3) and non-integrable 
(187r/31, 177r/97) triangles considered, the initial decay 
is exponential while the asymptotic decay is a power law, 
r(n) ~ n~^, with /3 = 1.035 in the integrable case and 
/3 — 1.085 in the non-integrable (NI) case. Thus pre- 
exponential decays are not exceptional and can persist 
for a long time in marginally stable systems. 

A more realistic situation should however include 
noise. For instance, imperfections can give rise to maps 
of the type considered in eq. (|l|). This leads to inter- 
esting results. For Gaussian noise (in Birkhoff momen- 
tum) with a = 0.000001, a single exponential decay dom- 
inates the survival fraction in the NI case for n > 5000 
while in the integrable case, the transition continues be- 
yond n — 14000. Moreover, the closely spaced eigen- 
values leads to a quasi-algebraic decay in the interval 
6000 < n < 14000 for the integrable case. Note that most 
trajectories remain largely unaffected for several hundred 
bounces for the value of a considered so that r(n) closely 
follows the noiseless case initially. Thus, even when the 
asymptotic decay is exponential, transition times can be 
very large. 

In the weak noise case, the evolution operator can be 
approximated as 



(14) 



where ^(x) = (/)(a;)/ |/'(/- ^(a;))|. Using a polynomial 
basis p] , a matrix representation of the operator can be 
constructed where the elements £„„ = ^^{£ ° (^)}- 
Writing a as ai+a^ where cto — 0^, a perturbation calcu- 
lation shows that the eigenvalues (and thus the difference 
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between Aq and Ai) decrease as cr^ when is small. The 
transition time therefore decreases with noise. 
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FIG. 3. r(n) (log scale) for the (7r/2, 7r/3, 7r/6) (inte- 
grable) and the (IStt/SI, 177r/97) (non-integrable) triangle 
along with the best fit exponential (dashed line) in the in- 
tegrable case. Here a = 0.00005 and Ago ~ 0.005. 

For a — 0.00005 (see fig. ^), the transition time decreases 
significantly and exponential decay sets in for n > 3500 
in the integrable case while the leading eigenvalue dom- 
inates from the beginning in the non-integrable case. 
Thus the gap between Aq and Ai increases with noise. 

In conclusion, the broad picture that emerges from 
these numerical experiments is as follows : 

(i) Additive noise makes the spectrum of the evolution 
operator discreet. 

(ii) When the dynamics is intermittent or regular and 
the noise weak, exponential decays may emerge only 
asymptotically due to the presence of closely spaced 
eigenvalues around the leading eigenvalue, Aq. These are 
remnants of the continuous spectrum that exists in the 
zero noise case. The transition phase in such a situation 
can mimic an algebraic decay. 

There are important fallouts of this conclusion. In 
experimental situations where noise is inevitable, signa- 
tures of exponential decays are not necessarily indicative 
of chaotic dynamics. For instance, semiclassical theory 
links Lorentzian line shapes observed in experiments on 
ballistic transport in chaotic microstructures to the ex- 
ponential decay in the survival fraction |jl^ . The present 
analysis however shows that noisy intermittent dynam- 
ics can also give rise to Lorentzian line shapes and it is 
interesting to note that there are instances where obser- 
vations on regular or marginally stable cavities have been 
found to be no different from chaotic cavities [0| . 
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s 


Ao 


A/it 


0.5 


0.99470 


0.99467 


0.7 


0.99508 


0.99502 


0.9 


0.99553 


0.99546 



TABLE I. The leading eigenvalue, Aq computed by discree- 
tizing C{x,y) compared to the best fit value A/it (see fig. 1) 
for three different parameter values. 
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